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O : ABSTRACT 

a 

We analyze V, I and H band HST images and two seasons of i?-band mon- 
i. itoring data for the gravitationally lensed quasar SDSS0924+0219. We clearly 

j_j ■ see that image D is a point-source image of the quasar at the center of its host 

^ galaxy. We can easily track the host galaxy of the quasar close to image D be- 

cause microlensing has provided a natural coronograph that suppresses the flux of 
the quasar image by roughly an order of magnitude. We observe low amplitude, 
uncorrelated variability between the four quasar images due to microlensing, but 
no correlated variations that could be used to measure a time delay. Monte Carlo 
models of the microlensing variability provide estimates of the mean stellar mass 
in the lens galaxy (O.O2M0 < (M) < l.OM©), the accretion disk size (the disk 
temperature is 5 x 10^ K at 3.0 x 10^^ cm ^ < 1.4 x 10^^ cm), and the black 
hole mass (2.0 x IO^Mq < Mbh Vo.i^ {L/LeY^^ < 3.3 x lO^M©), all at 68% con- 
fidence. The black hole mass estimate based on microlensing is consistent with 
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an estimate of Mbh = (7.3 ± 2.4) x lO'^'M© from the Mg II emission line width. If 
we extrapolate the best-fitting light curve models into the future, we expect the 
the flux of images A and B to remain relatively stable and images C and D to 
brighten. In particular, we estimate that image D has a roughly 12% probability 
of brightening by a factor of two during the next year and a 45% probability of 
brightening by an order of magnitude over the next decade. 

Subject headings: cosmology: gravitational lensing - microlensing - stellar masses 
- quasars: individual (SDSS0924+0219) - accretion disks - dark matter 



Introduction 



In the standard Cold Dark Matter (CDM) galaxy model, early-type galaxies are com- 
posite objects. Their mass is dominated by an extended dark matter halo tha t surrounds the 
lumin ous stars of the visible galaxy; any remaining baryons are left as hot gas (jWhite fc Rees 
19781 ). The halos grow by mergers with other halos, with a small fraction of the accreted 
smaller halos surviving as satellite halos ( "su b-structure" ) orbiting in the larger ha l os, but the 
mass fraction in these satellites is uncertain (IMoore et al.lll999l . iKlypin et al.lll999l . iGao et al. 



2004J . iTavlor fc Babulll2005 



and Zentner et al 



20051 ) . There is increasing evidence from time 



delay measurements (e.g. iKochanek et al.ll2005h and stellar dynamical observations (e .g. lo- 
cally, lR«nmi]wsky_eLalJl2003|, and in lenses iTreu fc Koopmansll2002l . iTreu et al.ll2006l ). that 
the density structure of some early-type galaxies on scales of l-2i?e is heterogeneous, but this 
needs to be changed from a qualitative assessment to something more quantitative. Thus the 
detailed balance between stars, dark matter and substructure (luminous or dark) remains a 
matter of debate. 

One approach to addressing these problems is to monitor variability in gravitational 
lenses. Variability in lenses arises from two sources. Correlated variability between the 



^ Based on observations obtained with the SmaU and Moderate Aperture Research Telescope System 
(SMARTS) 1.3m, which is operated by the SMARTS Consortium, the Apache Point Observatory 3.5-meter 
telescope, which is owned and operated by the Astrophysical Research Consortium, the WIYN Observatory 
which is owned and operated by the University of Wisconsin, Indiana University, Yale University and the 
National Optical Astronomy Observatories (NOAO), the 6.5m Magellan Baade telescope, which is a collab- 
oration between the observatories of the Carnegie Institution of Washington (OCIW), University of Arizona, 
Harvard University, University of Michigan, and Massachusetts Institute of Technology, and observations 
made with the NASA/ESA Hubble Space Telescope for program HST-GO-9744 of the Space Telescope Sci- 
ence Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. 
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images due to fluctuations in the source flux allows the measurement of the time delays At 
between the quasar images, which constrain the combination At oc {1 — {k))/Hq of the surface 
density near the lensed images (k) = (S)/Sc and the Hubble constant Hq to lowest order 
(Kochanek 2OO2I ) 1^ By measuring the surface density near the lensed images, we can strongly 
constrain the radial mass profile of the lens. Uncorrelated variability between the images 
is a signature of microlensing by stars in the lens galaxy. Microlensing can constrain the 
mass distributio n because the statistics of the va riations depend on the fraction of the local 
density in stars (ISchechter fc Wambsganssl l2002l ). In addition to providing an estimate of 
the surface density in stars near the lensed images, microlensing can also be used to estimate 
the mean stellar mass in the lens and to determine the structure of quasar accretion disks. 
Understanding microlensing is also required to improve estimates of the substructure mass 
fraction in radio-quiet lenses, because most other quasar source components are affected 
by both substructure and microlensing. In order to use time variability in lenses to probe 
these astrophysical problems, we have undertaken a program to monitor roughly 25 lenses in 
several optical and near-IR bands. The first results of the program and a general description 
of our procedures are presented in iKochanek et al.l (120051 ). 



In this paper we study the four-image Zs = 1.52 quasar lens SDSS 0924+0219 (llnada et al 



20031 ) using V, I and H-hand Hubble Space Telescope (HST) observations and the results 
from two se asons of monitorin g the system in the i?-band. The lens galaxy is a fairly isolated 
zi = 0.393 (lOfek et al.l l2005l ) early-type galaxy. The most remarkable feature of this lens 
is that it shows a spectacular flux ratio anomaly between the A and D images. These two 
images are merging at a fold caustic with a flux differenc e of nearly 3 magnit u des when they 



should have approximately equal fluxes by symmetry (IKeeton et al.l l2005l ) . iKeeton et al. 



(120061 ) recently measured the flux ratios in the Lya line and the adjacent continuum, finding 
that the anomaly is weaker in the emission line but still present. This indicates that the 
anomaly is partly due to microlensing, since the expected size difference between the broad 
line region and the optical continuum emission region should matter for microlensing but 
be irrelevant for substructure. Unfortunately, the continued existence of the anomaly in the 
emission line means either that the broad line region is not large enough to eliminate the 
effects of microlensing or that the flux ratio anomaly in SDSS0924+0219 is due the combined 
effects of microlensing and substructures. 

We present the HST data in ^ as well as a series of mass models for the system 
consisting of the observed stellar distribution embedded in a standard dark matter halo. In 



^The dimensionless surface density k is the surface density S divided by the critical surface density 
Ec = c^Doii/AnGDoi^Di^a,, where -Dql, £'os a-nd Z^ls a-re the angular diameter distances between the 
Observer, Lens and Source. 
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Owe present the light curves and model them using the Monte Carlo methods of iKochanek 
(|2004J ). The analysis allows us to estimate the mean stellar mass in the lens galaxy, the size 
of the quasar accretion disk and the mass of the black hole powering the quasar. In §3.41 
we present our predictions for the expected variability of this source over the next decade. 
Finally, we summarize our findings in §11 All calculations in this paper assume a flat ACDM 
cosmology with h = 0.7, fin = 0.3 and = 0.7. 



2. HST Observations and Mass Models 

In this section we discuss the HST observations and our mass model fits to the astro- 
metric and photometric measurements. 



2.1. HST Observations 



We observed the lens in the V- (F555W), /- (F814W) and i^-bands (F160W) using 
HST. The ~ 4380 sec \^-band and the ~ 4600 sec J-band images were obtained as eight 
dithered sub-images with the Wide Field Channel (WFC) of the Advanced Camera for 
Surveys (ACS) on 2003 November 18. The 5120 sec if-band image was obtained as eight 
dithered sub-images on 2003 November 23 using the Near- Infrared Camera and Multi-Object 
Spectrograph (NICMOS). The ACS data were reduced using the pyraf m ultidrizzle package, 
and the NICMOS data were reduced using nicred (see iLehar et al.ll2000l ). We focus on the 
results from the new i^-band image, shown in Fig. [H The four quasar images, the central 
lens galaxy and an Einstein ring image of the quasar host galaxy can easily be seen. 

We fit the if-band image using a photometric model consisting of four point sources 
and a de Vaucouleurs model for the lens galaxy. We modeled the lensed host galaxy as an 
exponential disk, using a singular isothermal ellipsoid (SIE) mass model for the lens galaxy. 
We also tried a de Vaucouleurs lens galaxy mass model, but doing s o caused no signi ficant 
change in the quality of the lensed host galaxy photometric fit (see iPeng et al.l |2006| for a 
discussion o f systematic issues in modeling lensed host galaxies). The fits were done with 
imfitfits (see lLehar et al.ll2000l ) using a range of bright PSF models, with the PSF producing 
the best overall fit being adopted for the final results. We then use the if-band fit as our 
reference, and we hold the astrometry and model structure fixed for the V- and J-band 
photometric fits. In Table (H we present astrometric and photometric measurements for 
the system. Our fi ts to the lens astrome t ry an d photometry are consistent with those of 
Keeton et all fl2006h and lEigenbrod et all J2005h . 
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The lens in SDSS 0924+0219 is an early-type galaxy with effective radius Re = 0'.'31 ± 
0'.'02, axis ratio q = 0.92 ± 0.02, major axis position angle 9^ = -27 ± 8°(East of North) 
and colors V — I = 1.44 ± 0.07 and V — H = 3.60 ± 0.05. Following the technique of 



Rusin fc Kochanekl (120051 ). we performed a Fundamental Plane (FP) analysis of the lens to 
determine if its observed properties fall on th e track of the exp ected Mass-to-Light (M/L) 
evolution with redshift for elliptical galaxies (IRusin et al.ll2003l ). Assuming that a nominal 
early-type galaxy undergoes a starburst phase at high redshift and then evolves passively 
thereafter, we expect elliptical galaxies to show a steady increase in M/L with decreasing 
redshift. If the SDSS0924+0219 lens galaxy is to lie on the present-day FP, then it requires 
a Alog(M/L) evolution of —0.42 ± 0.03 from z = 0.393 to z = 0. This is steeper than the 
mean value for lens galaxies around that redshift, Alog(M/L) = —0.22 ± - 04, im plying 
that the lens has a sr aaller than average M/L. Note, however, that iTreu et al.l (120061 ) find a 
faster evolution than iRusin et al.l (120031 ). which would be consistent with the value we find 
for SDSS0924+0219. 

One remarkable feature in this lens is how the stars (or a small satellite) in the lens galaxy 
have provided a natural coronograph at the location of image D. The quasar flux is suppressed 
by roughly an order of magnitude, making it very easy to see into the central regions of the 
host galaxy. We obtain a good fit with a host having a scale length of O'.'ll ± O'.'Ol, axis 
ratio of 0.74 it 0.05 and magnitude of H = 20.56 ± 0.14 mag. For comparison, we estimate 
an unmagnified magnitude of H = 20.40 ± 0.20 mag for the quasar. We also attempted 
a de Vaucouleurs model fit to the host galaxy, but we found that doing this resulted in a 
negligible change in the overall quality of fit. We are not able to discriminate between h ost 
galaxy models. We extracted the Einstein ring curve of the lens ( Kochanek et ahlboOll ) to 
use as one of the constraints on the mass models. 

Images A, B and C have similar I — H and V — I colors, while image D is significantly 
redder in I — H but of similar color m V — I. This is similar to the expected pattern 
for dust extinction, since an 0.5 mag difference m I ~ H should correspond to a 0.15 mag 
difference m V — I for an Ry = 3.1 extinction curve shifted to the redshift of the lens. 
On the other hand, the Einstein ring image of the host galaxy does not show color trends 
near the D image, which strongly argues against dust as the origin of the color differences. 
Moreover t he significant differe nces between the A/D flux ratio in Lya as compared to the 
continuum (IKeeton et al.ll2006l ). means that microlensing must be a significant contributor 
to the anomalous flux ratio. 
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2.2. Macro Models and Substructure 



We modeled the lens galaxy as the sum of a de Vaucouleu rs model with scale length 
= C'Sl embedded in an NFW (iNavarro. Frenk fc Whitdll996l ) model with a break radius 
Tc = ICy/O. The de Vaucouleurs and NFW models were ellipsoids constrained by the axis 
ratio and orientation of the lens galaxy in the if-band images, and we included an external 
shear to model any additional perturbations from the lens environment or along the line 
of sight. We constrained the mass model with the astro metry of the quasar i mages and 
the Einstein r ing curye der ived from the if-band images (IKochanek et al.l l200ll ) using the 
GRAVLENS (IKeetonl 120011 ) software package. We required the NFW and de Vaucouleurs 
components to be perfectly concentric, but we permitted the combined model to move within 
0^.^01 of th e meas ured galaxy center in order to optimize the fit. As described in detail by 
Kochanekl (120051 ) . it is not possible to determine the radial mass profile of the lens using 
these constraints, although it can be done using other constraints such as a time delay or 
stellar velocity dispersion measurement. Given this degeneracy, we generated a sequence of 
models parameterized by < Jm/l < 1, the fractional mass of the de Vaucouleurs component 
compared to a constant M/L model {Jm/l = !)• As expected, there is no significant difference 
in how well models with < Jm/l < 1 fit the constraints with the exception of pure dark 
matter models {/m/l ^ 0.1) that predict a detectable, fifth or odd quasar image near the 
center of the Einstein ring. The convergence, shear and for the range of Jm/l at each 
image location is presented in Table [2J 

Although much of the anomalous A/D flux ratio must be due to microlensing based 
on the smaller anomaly observed in the emission lines, we explored the extent to which 
the anomaly could be created by small satellites of the lens galaxy rather than by stars. 
For these tests we modeled the main lens as a singular isothermal ellipsoid (SIE) and then 
added a low-mass pseudo-Jaffe model as a perturber. We assigned the perturber Einstein 
radii of either O'.'Ol or O'.'OOS and tidal truncation radii of O'.'l and 0'.'06, respectively. Fig. [2] 
shows the goodness of fit, x'^/^dof, as a function of the perturber's position, where we 
fit both the astrometric constraints from the quasar image positions and the flux ratios 
from the if-band HST data. While the more massive satellite with an Einstein radius of 
O'.'Ol has difficulty adjusting the flux ratio without violating the astrometric constraints, the 
lower mass satellite can do so if properly positioned. Such a satellite would have a mass 
of ~ 10~^ that of the primary lens, roughly corresponding to the mass scale of globular 
clusters. While the emission line flux ratios largely rule this out as a complete explanation 
for SDSS 0924+0219, substructure could explain the continued existence of an anomaly in 
the emission line flux ratios. We include this calculation as an illustration that substructure 
can lead to anomalies as extreme observed here. 
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3. Microlensing 



3.1. Monitoring Data and Microlensing 



We have obtained somewhat more than two seasons of -R-b and monitoring data fo r SDSS 
0924+0219. Our analysis procedures are described in detail in lKochanek et al.l (120051 ) . so we 
provide only a brief summary here. We measure the flux of each quasar image relative to a 
sample of reference stars in each frame. We keep the relative positions of the components 
fixed, using the HST astrometry for the lensed components, and derive the PSF model and 
quasar flux for each epoch by simultaneously fitting the lens and the reference stars. The 
PSF is modeled by 3 nested, elliptical Gaussian components. The galaxy is included in the 
model at a constant flux which we determine by fitting all the data as a function of the 
galaxy flux and then adopting the galaxy flux that produces the best fit to the complete 
data set. We confirm that the lens galaxy flux is approximately constant at each epoch by 
examination of the residual galaxy flux after subtraction of the best-fit models, and we find 
no evidence for variability during our three seasons of monitoring. 

Most of our observations were ob tained at the queue -scheduled SMARTS 1.3m using the 
ANDICAM optical/infrared camera ( DePoy et al. 2003 )1^ Additional observations were ob- 
tained at the Wisconsin- Yale- Indiana (WIYN) observatory using the WIYN Tip-Tilt Module 
(WTTM) 3, the 2.4m telescope at the MD M Observatory using the MDM Eight-K|, Echelle 
and RETROCAM @ J Morgan et al.l |2005| ) imagers, the 3.5m APO telescop e using Spicam 
and the 6.5m Magellan Baade telescope using IMACS (IBigelow et al.lll999l ). Images taken 
under seeing conditions worse than 1'.'5 were discarded. We also added the photometry from 



Inada et al.l (120031 ) to extend our baseline to nearly four years for the microlensing calcula- 



tions. The i?-band light curves are displayed in Fig. [3], and the data are presented in Table 

El 

In Fig. El we also show the HST \^-band photometry scaled to the best-fit i?-band 
monitoring magnitude of Image A on the observation date. In the HST data, image D is 
~ 1 mag fainter relative to image A than our estimate on nearly the same date. After 
considerable experimentation, we concluded that our flux for image D may be contaminated 
by image A, although we found no correlation between the A/D flux ratio and the seeing. 



■^http: / /www. astronomy.ohio-state.edu/ANDICAM/ 
^http://www. wiyn.org/wttni/WTTM_manual. html 
'^http:/ /www. astro. columbia.edu/ arlin/MDM8K/ 
http://www.astronomy.ohio-state.edu/MDM/RETROCAM 
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Nonetheless, we include estimates of seeing at each observation epoch in Table [31 In the 
calculations that follow we will use both our image D light curve as observed and an image 
D light curve shifted 1 mag fainter to agree with the HST flux ratio. The shift had little 
effect on our results in § I3.2H3.4I 

In our mass models, the longest expected delay for the system is ~ 10 days for Hq = 
72 km s~^ Mpc~^. We see little evidence for correlated variability between the images on 
these time scales, so we cannot measure the time delays. For the present study, it seems safe 
to simply ignore the time delays. 



3.2. Microlensing Models 

For each of our macro models, parameterized by Jm/l, we generated 8 random realiza- 
tions of the expected microlensing magnification patterns for each image. We used patterns 
with an outer dimension of 20-Re where i?e = 5.7 x 10^^ (M/M©)^/^ cm is the Einstein radius 
for the average mass star. The map dimensions were 8192^, so we can model source sizes 
down to 3 X 10~^Re. The stars used to create the patterns were drawn from a Salpeter 
IMF with a dynamic range in mass of a factor of one hun dred. We modeled the acc retion 



disk of the quasar as a standard, face-on thin disk model (IShakura &: SunyaevI Il973l ) with 
a scale length of = /^^(M/Mq)^/^ where the microlensing behavior is determined by the 
source size scaled by the mean mass of the microlenses, fg. We have cho sen to use a thin 



disk m odel because it provides a context for interpreting the results, but iMortonson et al. 



( I2OO5I ) have shown that can simply be interpreted as the typical half-light radius for any 
choice of emission profile. We do, however, neglect the central hole in the disk emissivity to 
avoid the introd uction of an addit ional parameter. We fit the light curves using the Monte 



Carlo method of iKochanekl (120041 ). In this method, large numbers of trial light curves are 
randomly generated and fitted to the observed light curves. Bayesian statistical methods 
are used to combine the resulting distributions of values for the fits to the light curves to 
obtain probability distributions for the model parameters. 

We are interested in models where microlensing is responsible for any deviations of the 
image flux ratios from the lens model, so we assumed that the flux ratios of the macro model 
were correct up to a systematic uncertainty of 0.05 mag for images A-C and 0.1 mag for 
image D. We also allowed for an 0.02 mag systematic uncertainty in the photometry of images 
A-C and an 0.1 mag systematic uncertainty in the photometry of image D. These errors were 
added because the point-to-point scatter in the light curves is somewhat broader than the 
formal uncertainties in the photometry. With these assumptions, we have no difficulty finding 
light curves that fit the data well, with x^/^dof ~ 1-0. We generated 10^ trial light curves 
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for each source size, magnification pattern and mass model. Several example light curves 
which provide good fits to the data are shown in Fig. |H For each reasonable fit to the light 
curves, defined by x^/^dof < 2.3, we also generated an extrapolated light curve extending 
for an additional ten years beyond the last data point of the third season (2005 December 
14; HJD 2453719). We also repeated all the calculations shifting the image D light curve 
1 mag fainter to match the HST observations, finding few changes in the results. 

In order to convert the results from Einstein radius units, where all physical scales 
depend on the mean mass of the microlenses (M/Mq), we must assume either a probability 
distribution for the actual velocities or a prior for the mean stellar mass. Our velocity model 
includes the 176 km s^^ projected velocity of the CMB dipole onto the plane of the lens, a 
probability distribution for the one-dimensional peculiar velocity dispersion of galaxies at zi 
with rms value of 164 km s~^ and a one-dimensional stellar velocity dispersion in the lens 
galaxy of 219 km s~^ bas ed on the Einstein radius of the lens and an isothermal lens model. 



As discussed in detail by iKochanekl (120041 ). the lens galaxy pecuhar velocity dispersion and 
stellar velocity dispersion estimates are dependent upon the selected cosmology. We have 
chosen the standard flat concordance cosmology {Qm = 0.3 and Qa = 0.7) for our estimates. 
We also consider the consequences of using a limited range for the mean stellar mass of 
0.1 < (M/Mq) < 1.0. 

Figures [5] and [6] show the results for the system motions and the lens galaxy stellar mass 
estimate. The large flux ratio anomalies combined with the limited amount of observed 
variability means that the effective velocity in the system must be relatively low. We flnd 
that 280kms~^ ^ "^e ^ 749kms~^ (68% confldence), and this changes little if we adjust the 
D image light curve to be 1 mag fainter. If we compare the effective velocity distribution 
to our model for the possible distribution of physical velocities (Fig. [5]), we can estimate 
the mean microlens mass since the two velocities are related by Ve = Ve/ (M/Mq)^^"^. The 
broad range permitted for Ve translates into a broad range for the stellar mass estimates, 
with 0.02 M0 < (M) < 1.0 M0 (68% confldence). The low mass solutions correspond to 
large sources with high effective velocities and the high mass solutions correspond to small 
sources with low effective velocities. 

We were somewhat surprised to flnd that the present data do n ot distinguish between the 



lens models at all because we had based our expectations on the (ISchechter fc Wambsganss 



2OO2I ) picture in which models with low dominate the probability of flnding a faint saddle 
point image like D. We found instead that the probability distribution for /m/l is basically 
flat. This result is little affected by imposing the prior on the permitted mass range of the 
microlenses or by adjusting the image D light curve to be 1 r nag fa inter. In comparison, the 



microlensing models for SDSS 0924+0219 by iKeeton et al.l (120061 ) strongly favored models 
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with low ^^^,/^^. We believe the differences between the results are due to our use of finite- 
sized sources, which significantly enhance the probability of large demagnifications relative 
to large magnifications because the high magnification regions (caustics) are more affected 
by finite source sizes. Another source of differences is that we are analyzing a more strongly 
constrained problem by requiring that the models fit the observed light curves rather than 
simply fit the instantaneous flux ratios. 



3.3. Quasar Structure 



One objective of our monitoring program is to study the structure of quasar accre- 
tion disks. We start b y estimating the b lack h ole r nass using the Mrh, M k H line width, 
luminosity relations of iMcLure fc Jarvid (120021 ) and iKollmeier et al.l (120051) . We measured 
the Mg II(2800A) l ine width in spectra o btained by Ofek et al.l ( 2005 ) following the pro- 
cedures detailed in iKoUmeier eTaP J2005h . and we estimated the magnification-corrected 
continuum luminosity a t 3000A, La(3000A), using a power law fit to our HST data. For the 

jratio n we find a black hole mass of Mbh = (6.3±1.5) x 10"^ M©, 
2005h calibration we find Mbh = (7.3 ±2.4) x 10^ Mq. We adopt 



McLure fc Jarvis J 2002 ) call 



and for the (IKoUmeier et al 



the estimate based on the IKoUmeier et al.l (120051 ) calibration. Similarly, we estimate that the 
magnification-c orrected bolometri c luminosity of the quasar is L^d = (2.7±1.3) x 10^^ erg s~^ 
where we follow iKaspi et al.l (2000) in assuming that L^oi ~ 9 x ALa(5100A). This bolometric 
luminosity corresponds to an accretion rate M = (O.48±O.24)r7|^|M0 yr^^, where r] = 0.1?7o.i 
is the radiative efficiency of the accretion. Fig. [7] summarizes these "classical" constraints 
on the quasar. 

The new constraint we obtain from the microlensing observations is on the size of 
the quasar, which we can also e stimate using our accretio n disk model and the observed 
flux. A standard thin disk model (IShakura &: Sunyaevlll973l ) radiates as a black body with a 
temperature profile of T = T<.(i?/rs)~^/^, and the scale length we measure should correspond 
to the point in the disk where the temperature corresponds to the rest-frame wavelength of 
the filter band pass. For our i?-band data (2770A in the quasar rest frame), our scale length 
corresponds to the point where Ts{rs) — 5.2 x 10^ K. If the viscous energy release is radiated 
locally and we are well removed from the Schwarzschild radius, then the disk temperature 
and scale length are related to the black hole mass Mbh and accretion rate M by 

aT^ = ^GMsnM/STirl (1) 



(IShakura fc SunyaevI Il973l ). so a measurement of constrains the product MbhM. One 
means of estimating is to simply compute what it must be to produce the observed R- 
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band flux. Again assuming a standard, face-on tliin disk, tlie emission profile is 

I{R)oc[exp{{R/rsf')-lY' (2) 



(IShakura &: SunyaevI Il973l ). Assuming that the disk radiates locally as a blackbody, we 
integrate this emission profile over the physical extent of the disk to estimate its specific 
luminosity Ljy,em in the rest frame. Incorporating the system's geometry and correcting for 
redshift effects, we convert L^^em to Fx^obs, the specific flux in the observed frame. We then 
solve Fx^obs for to yield 

rs, , = 2.83 X 10^^^^ ( ^] ( ^y^' io-o.2iM,^,,-i9) ^-l 
V COS i \ th J J 

where Mx,obs is the observed magnitude, Dos/fu is the angular diameter distance to the 
quasar in units of the Hubble radius and i is the disk inclination angle, assumed to have 
an average value (i) = 60°. For SDSS 0924+0219, we find an unmagnified HST /-band 
magnitude I = 21.24 ± 0.25 mag, yielding the scale radius at the redshifted center of the 
HST /-band, = 6.3 ± 1.6 x lO^'' cm. Assuming the scaling of thin disk theory, 

we estimate an /?-band disk size of , = 5.0 ± 1.3 x 10^^ cm. We show the resulting 

^R,obs 

constraint on MbhM in Fig. [3 

Microlensing tests this theoretical calculation because the amplitude of the microlensing 
variability is controlled by the projected area of the source that smooths the magnification 
patterns. Fig. [S] shows our estimate of the scaled /?-band source size f^, which is related to 
the physical source size by = fs{M / Mq)^/"^ . The source must be quite compact relative to 
the Einstein radius of the typical microlens, with an exact bound that is presently difficult to 
determine because of the limited level of observed variability. We face two technical problems 
in extending Fig. [H] to smaller source sizes. The first problem is that our analysis code is 
presently limited to 8192^ magnification patterns, so when using an outer dimension large 
enough to produce a reasonable statistical representation of the magnification patterns it is 
difficult to resolve such small scales. The second problem is that even if we could resolve the 
smallest scales, we would find that the probability distribution flattens and becomes constant 
at small scale lengths. This occurs because the differences between small smoothing lengths 
are detectable only during caustic crossings - if our light curves do not extend to within a 
source size of a caustic, there is little effect from using a still smaller source size. Complete 
convergence at small scales will require a light curve with caustic crossings. 

Despite these problems, we can estimate the physical source size of the accretion disk 
reasonably well because there is a fairly strong degeneracy between the scaled source size f^, 
the scaled velocity Ve and the microlens mass scale (M/Mq) in the sense that more compact 
sources require smaller scaled velocities which implies a larger microlens mass scale for the 
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conversion to the physical source size (see lKochanekll2004l ) . Fig. [9] shows the estimates of the 
physical size both with and without the prior on the microlens masses. Reassuringly, the 
results depend only weakly on the prior. Nonetheless, will adopt the results with the mass 
prior: 3.0 x 10^^ cm ^ < 1.4 x 10^^ cm at 68% confidence. This is consistent with our 
earlier estimate from the continuum fiux but is a weaker constraint on MbhM, as shown in 

Fig.m 



3.4. Expectations for the Future Behavior of SDSS0924+0219 

For each light curve which passed a threshold of x^/^dof < 2.3, we generated a light 
curve for a period of ten years beyond our most recent observation. For each image, we then 
tracked the maximum change in the brightness in both the positive and negative directions 
on 1, 3 and 10 year time scales. The normalized, cumulative distributions of these maximum 
changes are shown in Fig. [TUl One of the more obvious predictions of Fig. [TU] is that images 
A and B are likely to remain constant while images C and D are likely to become brighter. 
One of the original motivations of this study was to estimate the time scale on which the 
fiux ratio anomaly would vanish as D moved out of a low-magnification region and became 
brighter. Here we find an approximately 12% chance that it will brighten by at least a factor 
of 2 in the course of the next year and a roughly 45% chance that it will brighten by more 
than an order of magnitude over the next decade. For the separate calculation in which we 
lowered the fiux of all points on the image D light curve by +1 magnitude, we predict a 
9% probabihty of image D brightening by a factor of two during the next year, and a 53% 
chance that image D will brighten by a factor of ten during the next ten years. 

Our expectation that D brightening would be the means of solving the anomaly was 
based on the preconception that D was a de-magnified sad dle point in a model with a smal l 



ratio between the stellar and total surface densities k^/k. ISchechter fc Wambsgansj (120021 ) 
demonstrated that in this regime there is an appreciable probability of strongly de- magnifying 
saddle point images like D. Our results confirm this finding. 



4. Conclusions 

During the course of our monitoring campaign we have observed uncorrelated variabil- 
ity in the four images of SDSS 0924-1-0219, evidence that microlensing is occurring in this 
system. Furthermore, our models demonstrate that microlensing is a viable explanation for 
the system's anomalous fiux ratios. This study does not rule out the alternative hypothesis 
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that dark matter substructure contributes to the anomaly, but it does firmly establish the 
presence of microlensing variability and the ability of microlensing to explain the anomaly. 
As we continue to monitor SDSS 0924+0219, we expect to eventually measure the time 
delay, thereby restricting the range of permissible halo models, and to steadily reduce the 
uncertainties in the estimated mean stellar mass, accretion disk structure and black hole 
mass. At some point over the next few years, we should also see dramatic changes in the 
fluxes of the merging images. 

We can also begin to con apare microlensin g estimates of the structure of quasar accretion 



disks. In our original study (iKochanekl |2004J ) . we modeled the significantly more luminous, 
but very similar redshift, quasar Q2237+0305 (My = -25.8 ± 0.5 versus My = -21.7 ± 0.7 
after correcting for magnification). As we would expect from accretion disk theory, the 
microlensing analyses indicate that the more luminous quasar has a significantly larger scale 
{ts — 4.1 X 10^^ cm versus ^ 6.9 x 10^^ cm) and black hole mass {Mbh — 1-1 x 10^ M© 
versus 1.3 x 10^ Mq). The next step is to combine the microlensing analyses of many lenses to 
explore these correlations in detail and to use the wavelength dependence of the microlensing 
variability to study the structure of individual disks. This next step should be possible very 
shortly. 
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Table 1. HST Astrometry and Photometry of SDSS0924+0219 



Component 


Astrometry 




Photometry 






ARA 


ADec 


H=F160W 


I=F814W 


V=F555W 


A 


= 


= 


17.96 ±0.02 


18.77 ±0.05 


19.61 ±0.01 


B 


+(y.'06i ± (y/003 


-l'/805 ± (y.'oos 


18.87 ±0.03 


19.64 ±0.07 


20.36 ± 0.05 


C 


-0'/968 ± 0'/003 


-0'/676 ± 0'.'005 


19.22 ±0.02 


20.22 ±0.09 


20.74 ±0.09 


D 


+o'.'536 ± c/ooa 


-0'.'427 ± 0"003 


20.64 ±0.06 


22.00 ±0.22 


22.94 ±0.13 


G 


-C'lSS ± 0'.'004 


-0'.'858 ± 0'/004 


17.23 ±0.04 


19.39 ±0.06 


20.83 ±0.03 



Table 2. Macroscopic Lens Mass Models 



/m/z. 




Convcrgi 


cncc K 






Shear 7 














A 


B 


C 


D 


A 


B 


C 


D 


A 


B 


C 


D 


0.1 


0.74 


0.71 


0.77 


0.75 


0.24 


0.21 


0.29 


0.28 


0.020 


0.018 


0.024 


0.021 


0.2 


0.67 


0.65 


0.70 


0.67 


0.30 


0.26 


0.37 


0.36 


0.042 


0.038 


0.050 


0.044 


0.3 


0.61 


0.59 


0.64 


0.61 


0.35 


0.30 


0.44 


0.43 


0.061 


0.055 


0.072 


0.062 


0.4 


0.53 


0.51 


0.56 


0.53 


0.42 


0.36 


0.53 


0.52 


0.11 


0.10 


0.13 


0.11 


0.5 


0.47 


0.45 


0.49 


0.46 


0.48 


0.41 


0.61 


0.60 


0.15 


0.14 


0.17 


0.15 


0.6 


0.40 


0.38 


0.44 


0.41 


0.54 


0.45 


0.69 


0.67 


0.19 


0.18 


0.23 


0.20 


0.7 


0.31 


0.29 


0.35 


0.31 


0.62 


0.51 


0.80 


0.77 


0.34 


0.32 


0.38 


0.34 


0.8 


0.25 


0.23 


0.29 


0.26 


0.67 


0.55 


0.88 


0.84 


0.45 


0.42 


0.50 


0.46 


0.9 


0.17 


0.16 


0.22 


0.18 


0.73 


0.61 


0.96 


0.92 


0.72 


0.70 


0.77 


0.73 


1.0 


0.13 


0.12 


0.17 


0.13 


0.77 


0.64 


1.02 


0.98 


1.00 


1.00 


1.00 


1.00 



Note. — Convergence k, shear 7 and the fraction of the total surface density composed of stars 
at each image location for the series of macroscopic mass models. 



Table 3. SDSS0924+0219 Light curves 



HJD X^/Ndof Seeing QSO A QSO B QSO C QSO D (Stars) Source 



2957, 


.814 





.70 


1^22 


2 


.812 it 0, 


.028 


3 


.515 it 0, 


.033 


3, 


.883 it 0, 


.054 


5 


.041 it 0, 


.217 


— 0, 


.032 it 


.003 


SMARTS 


2964, 


.789 


0, 


.65 


1'.'42 


2 


.864 ± 0, 


.024 


3 


.516 ± 0, 


.024 


3, 


.895 ± 


.042 


4, 


.840 ± 0, 


.144 


0, 


.004 ± 0, 


.003 


SMARTS 


2976 


.775 


0, 


.64 


1'.'31 


2 


.800 ± 0, 


.020 


3 


.464 ± 0, 


.022 


3, 


.892 it 


.039 


5 


.006 ± 0, 


.149 


— 0, 


.001 ± 


.003 


SMARTS 


2984, 


.721 


0, 


.56 


1'.'20 


2 


.811 ± 0, 


.036 


3 


.453 ± 0, 


,040 


3, 


.852 it 0, 


.068 


4, 


.632 ± 0, 


.189 


— 0, 


.037 ± 0, 


.003 


SMARTS 


2991 


.791 





.79 


1'.'08 


2 


.766 it 0. 


.016 


3, 


.531 it 0. 


.019 


3. 


.991 it 0, 


.033 


5 


.138 it 0. 


.128 


0. 


.012 it 0, 


.003 


SMARTS 


2995, 


.932 


1, 


.82 


1^20 


2, 


.785 it 0. 


.012 


3, 


.507 it 0. 


.013 


3. 


.928 it 0, 


.022 


5, 


.147 it 0. 


.105 


0. 


.084 it 0, 


.002 


MDM-EIGHTK 


2998, 


.787 


1, 


.31 


0{'94 


2, 


.804 it 0. 


.014 


3, 


.541 it 0. 


.017 


4. 


.047 ± 0, 


.029 


5, 


.258 it 0. 


.120 


0. 


.027 it 0, 


.003 


SMARTS 


3022 


.797 


1 


.82 


1^'09 


2 


.882 it 0, 


.019 


3 


.417 it 0, 


.018 


4, 


.161 it 0, 


.037 


4, 


.937 it 0, 


.116 


0, 


.015 it 


.003 


SMARTS 


3029 


.751 


2 


.13 


0'.'90 


2 


.865 ± 0, 


.014 


3 


.386 ± 0, 


.015 


4, 


.187 ± 


.029 


5 


.092 ± 0, 


.092 


0, 


.021 ± 


.003 


SMARTS 


3036 


.673 





.70 


0'.'99 


2 


.780 ± 0, 


.016 


3 


.441 ± 0, 


.019 


4, 


.216 ± 


.042 


5 


.209 ± 0, 


.144 


0, 


.004 ± 


.003 


OA T A Torino 

SMARTS 


3047, 


.794 





.94 


1'.'31 


2 


.749 ± 0, 


.024 


3 


.408 ± 0, 


.024 


4, 


.331 it 


.063 


5 


.033 ± 0, 


.193 


— 0, 


.013 ± 


.003 


OA T A T^rnO 

SMARTS 


3055, 


.689 


0, 


.68 




2 


.819 it 0. 


.009 


3, 


.555 it 0. 


.012 


4. 


.206 it 0, 


.020 


5, 


.353 it 0. 


.071 


— 0. 


.025 it 0, 


.002 


TV T A 771 T T A TS.T 

MAGrjLLAJN 


3056 


.744 





.80 


1^11 


2 


.757 it 0, 


.016 


3 


.428 it 0, 


.017 


4, 


.201 it 0, 


.036 


5 


.152 it 0, 


.133 


0, 


.016 it 


.003 


SMARTS 


3064 


.680 


1 


.30 


I'.'OB 


2 


.836 ± 


,018 


3 


.381 ± 


,018 


4, 


.187 ± 


.039 


4, 


.948 ± 0, 


.119 


0, 


.010 ± 


.003 


OA T A T^rno 

SMARTS 


3065 


.805 





.61 


1'.'28 


2 


.832 ± 0, 


.018 


3 


.389 ± 0, 


.017 


4, 


.179 ± 


.039 


4, 


.560 ± 0, 


.091 


0, 


.080 ± 


.002 


MDM-ECHELLE 


3071 


.649 


0, 


.53 


1'.'08 


2 


.800 ± 0, 


.033 


3 


.423 ± 0, 


.039 


4, 


.269 ± 0, 


.092 


4, 


.782 ± 0, 


.206 


— 0, 


.044 ± 0, 


.003 


SMARTS 


3075 


.774 


4, 


.24 


1'.'07 


2 


.828 it 0. 


.012 


3, 


.491 it 0. 


.013 


4. 


.005 it 0, 


.022 


5 


.046 it 0, 


.091 


0, 


.053 it 0, 


.002 


MDM-EIGHTK 


3078, 


.586 





.57 


I'.'li 


2, 


.791 it 0. 


.019 


3, 


.446 it 0. 


.020 


4. 


.234 ± 0, 


.045 


4, 


.888 it 0. 


.128 


0. 


.007 it 0, 


.003 


SMARTS 


3080 


.834 


0, 


.56 


1^'31 


2, 


.797 it 0. 


.017 


3, 


.419 it 0. 


.016 


4. 


.268 it 0, 


.041 


4, 


.742 it 0. 


.098 


— 0. 


.002 it 0, 


.003 


APO 


3087 


.648 


1 


.49 


0^'86 


2 


.770 it 0, 


.007 


3 


.408 it 0, 


.008 


4, 


.294 it 0, 


.012 


5 


.384 it 0, 


.043 


0, 


.188 it 0, 


.003 


WIYN-WTTM 


3089 


.599 


1 


.11 


1'.'04 


2 


.781 ± 0, 


.016 


3 


.409 ± 0, 


.017 


4, 


.198 ± 


.036 


5 


.150 ± 0, 


.126 


0, 


.011 ± 


.003 


OA T A T^rT^O 

SMARTS 


3090 


.666 


0, 


.79 


I'.'SO 


2 


.833 ± 0, 


.018 


3 


.336 ± 0, 


.013 


4, 


.450 ± 


.041 


4, 


.372 ± 0, 


.066 


0, 


.168 ± 


.003 


WIYN-WTTM 


3101 


.623 





.60 


1'.'20 


2 


.732 ± 0, 


.033 


3 


.417 ± 0, 


.037 


4, 


.420 ± 


.103 


4, 


.567 ± 0, 


.179 


— 0, 


.049 ± 


.003 


OA /t A T»nno 

SMARl S 


3108, 


.511 


2 


.18 


1'.'02 


2 


.710 it 0. 


.016 


3, 


.361 it 0. 


.018 


4. 


.213 it 0, 


.039 


5, 


.063 it 0, 


.129 


0, 


.009 it 0, 


.003 


SMARTS 


3115, 


.562 


1, 


.11 


1^19 


2, 


.697 it 0. 


.019 


3, 


.398 it 0. 


.020 


4. 


.273 ± 0, 


.049 


4, 


.889 it 0. 


.140 


0. 


.003 it 0, 


.003 


SMARTS 


3127, 


.473 


0, 


.88 


a'92 


2, 


.702 it 0. 


.020 


3, 


.364 it 0. 


.025 


4. 


.249 it 0, 


.058 


4, 


.735 it 0. 


.133 


— 0. 


.022 it 0, 


.003 


SMARTS 


3134 


.453 





.79 


0^'91 


2 


.631 it 0, 


.021 


3 


.467 it 0, 


.031 


4, 


.442 it 0, 


.080 


4, 


.861 it 0, 


.167 


— 0, 


.035 it 


.003 


SMARTb 




. OiU 


n 
U 




I'.'IO 


9 


-4- n 

.DOO ZH u. 


091 


Q 
O 


ziR7 -t- n 

.^D / III U. 


09^ 


A 
^, 


•^AS -t- fl 
.o^o III U 


.UDU 




.DOO ZH U. 


. ioi 


— u. 


.u±o ZH u 


. UUO 


oiVi/\r\.l o 


3159, 


.468 


0, 


.69 


1'.'07 


2, 


.596 ± 0, 


.024 


3, 


.484 it 0, 


.035 


4, 


.382 ± 0, 


.087 


5, 


.053 ±0, 


.238 


-0, 


.044 ± 0, 


.003 


SMARTS 


3169, 


.462 


0, 


.79 


1'.'28 


2, 


.635 ±0, 


.020 


3, 


.582 it 0, 


.025 


4, 


.269 ±0 


.057 


5, 


.135 ±0, 


.200 


-0, 


.012 ± 0, 


.003 


SMARTS 


3324, 


.813 





.42 


1'.'44 


2, 


.675 it 0. 


.030 


3, 


.882 it 0. 


.047 


4. 


.450 it 0, 


.099 


5, 


.528 it 0. 


.404 


-0. 


.032 it 0, 


.003 


SMARTS 


3327, 


.800 





.74 


1^'46 


2, 


.740 it 0. 


.027 


3 


.960 it 0. 


.042 


4. 


.429 it 0, 


.083 


5, 


.217 it 0. 


.265 


-0. 


.016 it 


.003 


SMARTS 


3330 


.766 





.74 


1('45 


2, 


.776 it 0, 


.037 


4, 


.050 it 0, 


.067 


4, 


.466 it 0, 


.118 


5, 


.129 it 0, 


.317 


-0, 


.038 it 0, 


.003 


SMARTS 


3344, 


.723 


0, 


.89 


1('43 


2, 


.767 it 0, 


.026 


3, 


.923 it 0, 


.039 


4, 


.628 it 0, 


.096 


5, 


.754 it 0, 


.399 


-0, 


.018 ± 0, 


.003 


SMARTS 


3349, 


.758 


0, 


.56 


1'.'49 


2, 


.703 ± 0, 


.032 


3 


.953 ±0, 


.048 


4, 


.971 ±0, 


.157 


5, 


.441 it 0, 


.384 


-0, 


.020 ± 0, 


.003 


SMARTS 


3357, 


.714 


1, 


.08 


1'.'26 


2, 


.765 it 0. 


.021 


3, 


.869 it 0. 


.029 


4. 


.586 it 0, 


.069 


5, 


.674 it 0, 


.297 


0, 


.003 it 0, 


.003 


SMARTS 


3360, 


.767 


1, 


.57 


I'.'U 


2, 


.799 it 0. 


.020 


3, 


.841 ± 0. 


.026 


4. 


.544 it 0, 


.056 


5, 


.306 ± 0. 


.192 


0. 


.007 it 0, 


.003 


SMARTS 


3376, 


.734 


1, 


.05 


1^15 


2, 


.780 it 0. 


.019 


3, 


.926 it 0. 


.029 


4. 


.533 it 0, 


.058 


5, 


.379 it 0. 


.202 


0. 


.006 it 0, 


.003 


SMARTS 


3382 


.701 





.61 


1^17 


2, 


.772 it 0. 


.022 


4, 


.012 it 0. 


.039 


4, 


.674 it 0, 


.082 


5, 


.468 it 0, 


.255 


-0, 


.017 it 


.003 


SMARTS 


3394, 


.741 


0, 


.61 


0'.'98 


2, 


.779 ± 0, 


.030 


3, 


.919 it 0, 


.064 


4, 


.908 it 0, 


.168 


5, 


.616 ±0, 


.426 


-0, 


.044 ± 0, 


.003 


SMARTS 


3403, 


.751 


1, 


.00 


0'.'98 


2, 


.801 ±0, 


.018 


3, 


.934 it 0, 


.028 


4, 


.733 it 0, 


.062 


5, 


.107 ±0, 


.141 


0, 


.003 ± 0, 


.003 


SMARTS 


3410, 


.781 


0, 


.99 


I'.'IS 


2, 


.767 ±0, 


.019 


3, 


.905 ±0, 


.028 


4, 


.681 it 0, 


.066 


5, 


.257 ±0, 


.183 


-0, 


.001 ± 0, 


.003 


SMARTS 


3427 


.737 





.55 


1'.'32 


2, 


.705 it 0. 


.041 


4, 


.025 it 0. 


.085 


5. 


.074 it 0, 


.247 


5, 


.181 it 0, 


.399 


-0, 


.031 it 0, 


.003 


SMARTS 


3435, 


.698 


1, 


.29 


(y.'99 


2, 


.793 it 0. 


.019 


3, 


.872 it 0. 


.029 


4. 


.666 it 0, 


.063 


5, 


.155 it 0. 


.158 


-0. 


.013 it 


.003 


SMARTS 


3443, 


.698 


1, 


.45 


1'.'29 


2, 


.783 it 0. 


.024 


3 


.900 it 0. 


.033 


4. 


.692 it 0, 


.080 


4, 


.858 it 0. 


.161 


-0. 


.006 it 0, 


.003 


SMARTS 
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Table 3 — Continued 



HJD 




Seeing 




QSO A 




QSO B 




QSO C 




QSO D 




(Stars) 




Source 


3450.662 


0.66 


1'.'20 


2, 


.701 ± 0.038 


3, 


.763 ± 0.066 


5, 


.006 ± 0.224 


6, 


.019 ± 0.807 


-0, 


,031 ± 0, 


.003 


SMARTS 


3461.643 


0.62 


1^46 


2 


.703 ± 0.029 


3 


.953 ± 0.045 


4, 


.965 ±0.141 


5, 


.626 ± 0.424 


-0, 


,025 ± 


.003 


SMARTS 


3468.550 


0.95 


1'.'28 


2 


.766 ±0.024 


3 


.930 ± 0.035 


4, 


,732 ±0.087 


5 


,106 ±0.200 


-0, 


,020 ± 


,003 


SMARTS 


3477.514 


0.49 


I'.'OB 


2 


.676 ±0.024 


3, 


.846 ± 0.044 


4, 


,604 ±0.099 


5 


,900 ± 0.482 


-0, 


,033 ± 


,003 


SMARTS 


3490.481 


0.77 


I'.'oe 


2, 


.670 ±0.019 


3, 


.825 ± 0.030 


4, 


,689 ±0.073 


4, 


,799 ±0.136 


-0, 


,010 ± 


,003 


SMARTS 


3502.497 


1.13 


1'.'06 


2, 


.632 ± 0.020 


3, 


.776 ± 0.031 


4, 


,662 ± 0.075 


4, 


.938 ±0.161 


-0, 


,013 ±0, 


.003 


SMARTS 


3520.515 


0.85 


1^34 


2, 


.520 ± 0.025 


3, 


.763 ± 0.038 


4. 


,739 ±0.112 


5, 


.261 ± 0.306 


-0, 


,024 ± 0, 


.003 


SMARTS 


3528.450 


0.77 


1^13 


2, 


.587 ± 0.027 


3, 


.743 ± 0.048 


4, 


,901 ±0.152 


5, 


.322 ± 0.343 


-0, 


,034 ± 0, 


.003 


SMARTS 


3676.836 


0.50 


1^40 


2, 


.185 ±0.020 


3 


.529 ± 0.033 


4, 


,584 ±0.114 


5, 


.174 ± 0.306 


-0, 


,030 ± 


.003 


SMARTS 


3686.836 


0.77 


1'.'22 


2, 


.120 ±0.015 


3, 


.518 ±0.025 


4, 


,375 ± 0.066 


5, 


.295 ± 0.258 


-0, 


,014 ±0, 


,003 


SMARTS 


3696.825 


0.54 


1'.'31 


2, 


.101 ±0.025 


3, 


.614 ±0.059 


4, 


,354 ±0.140 


5, 


,088 ± 0.403 


-0, 


,043 ± 0, 


,003 


SMARTS 


3704.817 


0.76 


I'.'ia 


2, 


.184 ±0.014 


3, 


.520 ±0.022 


4, 


,364 ± 0.057 


4, 


,969 ±0.164 


-0, 


,006 ± 0, 


,003 


SMARTS 


3711.988 


1.35 


0'.'99 


2, 


.221 ± 0.007 


3, 


.538 ± 0.008 


3, 


,909 ±0.011 


5, 


.206 ± 0.046 


0, 


,180 ±0, 


.003 


MDM-RETROCAM 


3712.791 


1.23 


1^03 


2, 


.213 ±0.013 


3, 


.505 ± 0.020 


4. 


,220 ± 0.044 


5, 


.099 ±0.162 


-0, 


,001 ± 0, 


.003 


SMARTS 


3719.804 


0.54 


1^13 


2, 


.173 ±0.019 


3, 


.586 ± 0.041 


4, 


,506 ±0.109 


5, 


.092 ± 0.279 


-0, 


,041 ± 0, 


.003 


SMARTS 



Note. — HJD is the Heliocentric Julian Day — 2450000 days. The goodness of fit of the image, /^doft is used to rescale the formal 
uncertainties when greater than unity (see text). The QSO A-D columns give the magnitudes of the quasar images relative to the comparison 
stars. The (Stars) column gives the mean magnitude of the standard stars for that epoch relative to their mean for all epochs. 
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B 



NICMOSHband(F160W) Residuals following subtraction 

Fig. 1. — if-band images of SDSS0924+0219. The left panel shows the original image. Note 
that image D shows the Airy ring of a point source and is markedly fainter than image A. 
The right panel shows //-band residual near each quasar image after subtracting the four 
quasar images and the lens galaxy to make the Einstein ring image of the quasar host galaxy 
more visible. The noisy residuals near each image are largely due to small errors in the PSF 
models. 
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Fig. 2.— Substructure models for the SDSS0924+0219 flux ratio anomaly. The left (right) 
panels show the regions where a pseudo-Jaffe model with an Einstein radius of (X'Ol ([/.'OOS) 
can remove the A/D flux ratio anomaly without signiflcantly worsening the constraints on 
the quasar image positions or the Einstein ring. The main lens is modeled as a singular 
isothermal ellipsoid (SIE) and the pseudo-Jafi^e models are tidally truncated at C'l ((/.'OG). 
The position of image D is indicated with a white point, contours are indicated in the 
flgure legends. 
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Fig. 3. — SDSS0924+0219 i?-band light curves for images A-D. The data points for image 
D are offset to improve their visibihty. The curves are derived from a joint, high-order 
polynomial fit for the source light curve combined with lower order polynomials for the 
microlensing variability of each image (see lKochanekll2005l for details). Symbols: Image A- 
solid squares, Image B-open triangles. Image C-open squares and Image D-solid triangles. 
The large solid triangles at 2962 days are the HST ^-band photometry referenced to Image 
A. 
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Fig. 4. — The five best fitting i?-band light curves extrapolated for an additional 10 years 
across their magnification patterns. In order to show the "goodness of fit" of the theoretical 
microlensing lightcurves, we plot as points the observed image flux minus our model for the 
intrinsic variability of the source. In most of these light curves we see a brightening of images 
D and C. 
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source plane (h ^<M/MQ>^/^km/s) 

Fig. 5. — Normalized probability distribution for the effective source plane velocity {ve, 
heavy solid line) as compared to our estimated probability distribution for the true source 
plane effective velocity (light solid line). Since the average microlens mass (M) is related 

1 /2 

to the two velocities by Vg oc Vg/ (M) ' , high (low) ratios of Ve/vg correspond to low (high) 
mass microlenses. 
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Fig. 6. — Probability distribution for the average stellar mass (M) in the lens galaxy. The 
uncertainty is relatively large because (M) oc but it is consistent with normal stellar 
populations. 
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Fig. 7. — Constraints on the mass and accretion rate of the SDSS0924+0219 quasar. The 
dashed horizontal hnes show the estimated black hole mass based on the Mg II 2800A emis- 
sion line width. The vertical dash-dotted lines show the estimated accretion rate assuming a 
radiative efficiency of 10% and the kasni et, ali H) method for estimating the bolometric 
luminosity. The dotted lines show the constraint on the product MbhM from the continuum 
flux and thin accretion disk theory. The shaded regions show the constraint on the product 
MbhM from the microlensing data with a prior of 0.1 < (M) < 1.0 on the mean 
mass of the microlenses. The accretion is super-Eddington in the lower right corner, so we 
terminated the constraints on the line where L = LEdd- AH values plotted are scaled to 
h = 0.7. 
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(h-i<M/MQ>i/2 cm) 

Fig. 8. — Probability distributions for the scaled source size fg. The heavy dashed line shows 
the estimate for fg including a prior of 0.1 Mq < (M) < I.QMq on the mass of the stars. 
The vertical line shows the Einstein Radius {Oe) of the average mass star. 
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Fig. 9. — Probability distributions for the physical source size r^. The dashed curve shows the 
estimate for with a prior of 0.1 Mq < (M) < 1.0 Mq on the mean mass of the microienses. 
The vertical hne shows the Schwarzschild radius Rbh — ^GMbh/c^ of a IO^Mq black hole. 
The last stable orbit for a Schwarzschild black hole is at 3 Rbh- 
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Am (mags) 

Fig. 10. — Fraction of light curves with 'X^ /Noof < 2.3 that undergo a change in magni- 
fication larger than Am magnitudes towards either brighter (right) or fainter (left) fiuxes 
after one (dotted), three (dashed) or ten (solid) years. The vertical lines show the largest 
observed Am for our present light curves. 



